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Abstract 

The anisotropy of the magnetic field dependence of the specific 
heat of YBa2Cu307_,5 can be used to identify different low-energy 
excitations, which include phonons, spin-i particles, and electronic 
contributions. With a magnetic field H applied perpendicular to the 
copper oxide planes, we find that the specific heat includes a linear- 
T term proportional to y/H. The nonlinear field dependence of the 
density of states at the Fermi level suggests that there are quasiparticle 
excitations throughout the entire vortex, not just in the vortex core. 
The ^/HT term agrees quantitatively with G. Volovik's prediction for 
a superconductor with lines of nodes in the gap. A similar, but much 
smaller, effect is predicted for fields parallel to the planes, and sensitive 
measurements of the in-plane anisotropic magnetic field dependence 
of the specific heat could be used to map out the nodes. 
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1 Introduction 



Among recent experiments supporting the possibility of lines of nodes in the 
gap function of some of the high-T c materials [l], measurements on single- 
crystal YBa2Cu307_<5 (YBCO) show a linear temperature dependence of the 
penetration depth 0, as predicted for superconductors with lines of nodes in 
the clean limit ||. The density of states N(E) in such superconductors rises 
linearly with energy at the Fermi level in zero field, which should result in a 
quadratic term aT 2 in the specific heat ||. A material with a finite N(Ep), 
such as a normal metal or a disordered superconductor with lines of nodes 
0, has a linear term jT in the zero-field specific heat, with 7 oc N(Ep). 

G. Volovik predicted that for a superconducting state with lines of nodes, 
the dominant contribution to the magnetic field dependence of the density 
of states at the Fermi level is N(Ep,H) = KN n JH/H c2 , where N n is the 
normal-state DOS and k is a factor of order 1 |J. The factor k is defined by 
the vortex lattice structure, by the slope of the gap near the gap node, and by 
possible structure in N n . The quasiparticles which contribute to N(Ep, H) 
are outside the vortex cores, and are close to the nodes in momentum space. 
By contrast, it is usually thought that in a clean s-wave superconductor the 
quasiparticles are confined to the vortex cores 0. 

In previous papers || [J we reported measurements of the magnetic field 
dependence of the specific heat of twinned and untwinned single crystals 
of YBa2Cu30"6.95 which were prepared identically to those showing linear 
temperature dependence of the penetration depth || [Uj . These crystals have 
a linear term j(0)T in the zero-field specific heat, a quadratic term aT 2 which 
is only a few percent of the total zero- field specific heat, and an additional 
field- dependent linear-T term 7j_(i?)T which obeys / y±(H)T = A^/HT. Two 
of these terms, aT 2 and A\^HT, agree quantitatively with the predictions of 
lines of nodes. 

In this paper we summarize the previous analysis, including alternative 
fits, and present results on the sample dependence of the zero-field linear 
term, j(0)T. We also discuss the physical origin of the dependence of 
the density of states and the prospects for observing this effect in a magnetic 
field applied parallel to the plane. In principle, measurements of the in-plane 
anisotropic magnetic field dependence of the specific heat could be used to 
map out the nodes. 
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2 Technique 



The specific heat was measured using a relaxation method described in de- 



tail elsewhere ||1 1|| . The sample was mounted on a sapphire substrate with 
a weak thermal link (thermal conductance k w ) to a constant-temperature 
copper block. At each temperature and field, k w was measured by apply- 
ing power P = K w (T samp i e — Tuock)- A smaller temperature difference was 
then used to make a relaxation measurement of the thermal time constant 
t = Ck w , where C is the total heat capacity. The precision of the mea- 
surement in the temperature range 2-10 K is ~ 0.5 %. The addenda heat 
capacity can be as high as half of the total heat capacity, limiting the abso- 
lute accuracy in some measurements to ~ 10 %. The addenda specific heat 
is mostly due to phonons, so that most of the systematic uncertainty effects 
only the Debye (3T 3 term. Measurements of an empty substrate showed that 
the addenda heat capacity does not depend on the applied magnetic field 
within the precision of the data. The accuracy of the field dependence of the 
specific heat is better than 1 %. 



3 Crystals 

All of the samples are single crystals of YE^C^C^.^ grown by the same flux- 
growth technique as crystals which show linear temperature dependence 



of the penetration depth ||. Results have been obtained on two nominally 
pure, twinned single crystals (samples Tl and T2). A 0.8 mg piece was 
cleaved from one of the twinned crystals, sample Tl, and detwinned to create 
sample Ul. About 20 % of the surface of sample Ul, near two of the corners, 
had a high density of visible twins remaining after the detwinning process. 
Sample T2 was mostly detwinned to create sample U2, which had visible 
twin boundaries remaining. Both Ul and U2 were reoxygenated subsequent 
to detwinning. 



4 Results 

Figure [I] shows the specific heat of sample Tl in an 8 Tesla field applied ei- 
ther parallel or perpendicular to the copper oxide planes. The total specific 
heat includes the Debye phonon contribution, f3T 3 , which has a similar value 
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in all three data sets; the zero-field linear term 7(0)T, which appears as a fi- 
nite intercept on the c/T-axis; a magnetic contribution from spin-- particles, 
ncschottk y (T, H), which appears as a low-temperature bump peaking around 
3-4 Kelvin in both parallel and perpendicular fields; and a substantially in- 
creased linear term in perpendicular field, r y±(H)T. 



4.1 Perpendicular Field — Global Fit 

To model the perpendicular field dependence, identical analyses are per- 
formed on data from samples Tl (H = 0, 0.5, 2, 4, 6, and 8 Tesla) and Ul 
(H = 0, 0.5, 2, 4, 6, 8, and 10 Tesla). A global fit is used to constrain the 
concentration of spin-| particles, n, and the phonon specific heat, /3T 3 , to 
be field-independent. Allowing for the possibility of a small quadratic term, 
aT 2 , as predicted for lines of nodes, the zero- field specific heat is described 
by: 

C (T, 0) = 7 (0)T + aT 2 + pT 3 . (1) 



The specific heat for H > 0.5 T is described by: 



(gfi B H/k B T), 



(2) 



where 

CSchottky(x) = X — 2 . (3) 

All data sets are fit simultaneously to Equations [1], and |^ with a Lande 
g- factor of g — 2 ||. Thus, (3 and n are constrained by multiple data sets, 
7(0) and a are constrained by the zero-field data, and j±(H)T is allowed a 
different value at each field. The statistical confidence levels are determined 
by a bootstrapping technique |I2" |, in which alternative data sets are ran- 



domly chosen from the original data set and the fitting procedure is repeated. 
This procedure for determining the statistical error bars avoids making any 
assumption about the statistical error distribution of the data set. 

The parameters determined from this fit are shown in Table [I], with 90 % 
statistical confidence levels. The coefficient of the Schottky anomaly, n, 
indicates 0.1% spin-| particles per Copper. The coefficient of the Debye 
term, f3 = 0.39 mJ/molK 4 , corresponds to a Debye temperature of 400 K. 
The difference in the Debye term between the two samples is within the 
systematic error of the Debye term of the addenda. The rms deviation of 
the data from the fit is 0.8 % for Tl and 2.7 % for Ul. The higher scatter in 
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the data for the untwinned crystal Ul results from the smaller sample size. 
The data and fits for both samples are shown in Figure |2], with the phonon 
specific heat (/3T 3 ) subtracted. 

Table 1: Global Fit Parameters (mJ, mol, K, and T) 





Sample Tl 


Sample Ul 


7(0) 


3.0 ±0.1 


2.1 ±0.1 -0.2 


n 


24 ± 1 


23 ± 1 


P 


0.392 ±0.001 


0.380 ±0.004 


a 


0.11 ±0.02 


0.10 ±0.06 


1±{H) 


0.91^ 





In the measurements on single crystals of YBCO reported here, the 
phonon specific heat obeys the Debye /3T 3 law up to 8 K. In the zero-field 
data set on sample Tl, for example, fits from 2 to 4 K and from 4 to 7 K 
give the same value of (3 within statistical error, indicating that higher pow- 
ers of temperature are not necessary to describe the data. Above 8 K, the 
deviations from /3T 3 are better described by a gapped excitation such as an 
optical phonon mode than by a sum of T 5 , T 7 , and T 9 terms. In contrast, 
the data on many polycrystalline samples require a large T 5 term to describe 
the phonon specific heat even below 5 K fi~3| . In order to avoid the extra 
parameters necessary to describe the phonon specific heat above 8 K, all of 
the fits reported in this paper are restricted to T < 7 K. 

The field-dependent linear term, 7_|_ (i?)T, has a nonlinear dependence 
on field (Figure |3|). This nonlinear field dependence is well-described by 
7±(#) = AVH, with A = 0.9 mJ/molK^T 1 / 2 for both samples Tl and Ul. 

4.2 Perpendicular Field — Alternative Fits 

The qualitative nonlinearity of J±(H) is robust to the assumptions used to 
describe the total specific heat, and is well described by J±(H) = A\^H. 
To check the assumptions of the global fit and the interdependence of the 
aT 2 and A\/~HT terms, each data set of sample Tl is fit independently to 
Equation |^. There are thus three parameters each for five data sets (H = 0.5- 
8 T) and two parameters for the zero-field data set (note that cs c hottky{T, 0) = 
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0), totaling seventeen parameters for six data sets. No zero-field quadratic 
term was allowed in this fit. A small quadratic term, if present, could easily 
be absorbed into a slightly larger cubic term, and this fit results in a (3{H = 0) 
which is ~ 5 % higher than (3{H ^ 0). The coefficients of the linear term are 
shown in Figure §a with error bars at the statistical 90 % confidence level. 
Fitting the coefficients of the linear term to j±(H) = Ay/HT gives A = 0.89 
mJ/molK^T 1 / 2 . 

In another check which avoids the zero-field data entirely, the five data 
sets with H > 0.5 T are globally fit to Equation ||, again keeping [3 and 
n independent of field. There are thus two global parameters, plus values 
of (7(0) + 7j_ (//")) in five fields (Figure $jp). The resulting linear-T term is 
better described by a y/~H dependence than by an H dependence, and from 
the fit shown in Figure |4]b, A = 0.87 mJ/molK 2 T 1//2 . This fit also returns an 
extrapolated value of 7(0) = 3.1 mJ/molK 2 , which is in agreement with the 
value determined by independently fitting just the zero- field data. 



4.3 Parallel Field 

The data on sample Tl in an 8 Tesla parallel field show a linear-T term, 
j\\(H)T, which is increased by about 0.5 mJ/molK 2 from 7(0), as determined 
by fitting the data in Figure |I] to Equation Q Because of uncertainties in the 
Schottky anomaly and the sample alignment, the results are also consistent 
with no increase in the parallel-field linear term, j\\(H)T. 



4.4 Zero-field Linear Term 

The coefficient of the zero-field linear term, 7(0)T, in the specific heat of 
each of the four samples is shown in Table 0. The zero-field linear term 
for the two twinned samples, samples Tl and T2, is 7(0) = 3 mJ/molK 2 , 
while in both untwinned samples, Ul and U2, the linear term is reduced 
to 7(0) = 2 mJ/molK 2 . In contrast, most YBCO samples have 7(0) > 4 
mJ/molK 2 [jnj. The coefficient of the linear term in the normal state of 
optimally-doped YBCO, j n , has been shown by other measurements to be 
about 20 mJ/molK 2 [ITS], ITJ]. The residual density of states N(Ep) in the 



superconducting state is therefore about 15 % of the normal-state density of 
states for the twinned samples and about 10 % for the untwinned samples. 



6 



Table 2: Zero-Field Linear Term, 7(0) (mJ/molK 2 ) 



Sample 7-5 Untwinned 7(0) N(E F )/N n 



Tl 


6.95 


no 


3.1 


~ 15% 


Ul 


6.95 


yes 


2.1 


~ 10% 




6.97 


yes 


1.9 


~ 10% 


T2 


6.95 


no 


2.8 


~ 15% 


U2 


6.95 


yes 


2 


~ 10% 



5 Discussion 

5.1 Comparison with Lines of Nodes 

The quadratic term is predicted to disappear in a magnetic field, where the 
energy dependence of the density of states close to the nodes, N(E, H = 
0) oc \E — E F \, is replaced by a finite N(E F , H). The aT 2 term is compa- 
rable in magnitude to the systematic uncertainty in the addenda T 3 phonon 
specific heat, but is several times larger than the relatively small systematic 
uncertainty of the field dependence, < 1 %. The global fit is thus the only 
possible way to identify the quadratic term in these data sets. The slight 
positive slope of the zero-field data set in Figure || can be described either 
as a phonon term which is about 5 % larger in zero field, which is difficult to 
explain, or as an aT 2 term. 

From the slope of the density of states in a superconductor with lines of 
nodes it is possible to predict a ~ 7„/T c « 0.2mJ/molK 3 for YBCO within 
factors of order unity ||, in good agreement with the value a = 0.1 mJ/molK 3 
obtained from the global fits (Table on samples Ul and Tl. 

For the field dependence of the density of states, Volovik predicts N(E F , H) = 
nN n ^J H / H c2 0, where k is of order 1 and is defined by the vortex lattice 
structure, by the slope of the gap near the gap node, and by possible struc- 
ture in N n . Taking y n = 20mJ/molK 2 |3j, |TJ and H c2 ^ = 150 T gives the 
prediction y±(H) = AyH with A = k * 1.6 mJ/molK 2 T 1//2 , in good agree- 
ment with the value A = 0.9 mJ/molK 2 T 1//2 found for both samples Ul and 
Tl. 

Both of the above comparisons used the measured value of j n , the angular 
averaged normal state density of states. Because both the aT 2 and Ay/HT 
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terms are a consequence of the nodes, these terms should be sensitive to the 
normal-state density of states close to the nodes. The above comparisons 
may indicate that the density of states close to the nodes is smaller than 
the angular average of the density of states, as indicated by photoemission 
measurements on BSCCO fL5|. 

These predictions do not distinguish between different types of lines of 



nodes, such as d-wave or extended s-wave [jlgj . In principle, this interpre- 
tation could also apply to a gap function with no nodes, but with a very 
small minimum gap, A <C ksT m 0.5 meV. In order to produce these re- 
sults, such a gap function would need to have an unusual energy depen- 
dence of the density of states similar to that associated with lines of nodes, 
N(E,H = 0) oc \E-E F \. 

A quadratic T 2 term has been known for some time in the specific heat 
of heavy-fermion superconductors, and has recently been reported in the 
specific heat of La 2 - a; Sr :r Cu04 [|17|. Two other recent works support the 
existence of a y/~H term in superconductors with lines of nodes: a reanalysis 
of existing data on polycrystalline YBCO samples shows a nonlinear j(H)T 
term which appears consistent with \fH fll8f , with similar but less dramatic 



curvature on Ca- and Sr-doped L^CuC^ samples, and a v H term has also 



been recently reported in the low-temperature specific heat of UPt3 [IS 



5.2 Residual Density of States 

The total residual density of states N(E F ) in the superconducting state, as 
determined from the 7(0)T term in the specific heat, is about 15% of the 
normal-state density of states for the twinned samples and about 10 % for 
the untwinned samples (Table 0). The residual density of states can also be 
extracted from fitting the temperature dependence of the penetration depth 
Ifjl to the d-wave model with scattering of Hirschfeld et al. ||. In samples 
similar to the samples measured here, these fits are inconsistent with a density 
of states which is greater than 1 % of N n ||. It appears difficult to reconcile 
the zero-field linear term 7(0)T with the aT 2 and A\fRT terms, and with 
the interpretation of the temperature dependence of the penetration depth. 
It is possible that the zero-field linear term has a separate origin from the 
excitations which give rise to the quadratic term and the field dependence of 
the specific heat. This assumption leads to a self-consistent analysis of the 
specific heat, and is given additional support by the measurements on the 
untwinned crystal Ul, in which 7(0) is decreased substantially but the aT 2 
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and Ay/~HT terms are unchanged within statistical error bars. 



5.3 Interpretation of the A\J HT term 

Conventional high-K superconductors are expected to have a linear-T term 
in the specific heat which is proportional to the volume of normal material in 
the cores, c ~ ■y n TH/H C 2 0, [HJ, and such an HT term has been observed in 
A15 superconductors ||21|| . The analysis which predicts this HT term does not 



apply to cuprate superconductors for two reasons. First, the small coherence 
length of the cuprates may result in a large excitation gap for quasiparticles 
which are confined to the vortex cores 0, p2| . Secondly, Caroli, deGennes, 
and Matricon |7| assumed a fully gapped superconductor in showing that the 
quasiparticle excitations are confined to the vortex core. 

Some form of a linear-H term might be expected in the specific heat 
for any excitations which are confined to the vortex cores. Although there 
may in principle be some dependence of the core excitations on magnetic 
field, it is likely that for H <C H c2 the core excitations would be roughly 
independent of field. Thus, any signal resulting from these excitations would 
simply be proportional to the number of cores, or linear in H. The nonlinear 
field dependence of the observed N(Ep, H) suggests excitations outside the 
vortex core. 

The physical origin of the density of states predicted by Volovik || is the 
same as the physical origin of the depairing current in a conventional super- 
conductor. In the presence of a supercurrent with velocity v S) the quasipar- 
ticle excitation spectrum is shifted by an amount k ■ v s . Far from the 
vortex core in a fully gapped superconductor, this shift is not large enough 
to change the density of states at the Fermi level, which is given by 

N(Ef) = / / - * • ■ ( 4 ) 

For a superconductor with lines of nodes, the shift is significant everywhere 
that the superfluid velocity is not zero. Assuming a vortex superfluid flow 
v s oc 1/r, the local density of states at the Fermi level also falls off as 1/r 
(Figure ^). Integrating over an entire vortex, with the intervortex spacing 
R{H) as the upper cutoff in the integral, gives a density of states per vortex 
N(Ep,H) oc R(H) oc 1/yH. Multiplying by the number of vortices, which 
is proportional to H, then gives the total density of states N(Ep, H) oc \[H. 
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The shift in the excitation spectrum depends on the angle between the 
local superfTuid velocity and the nodes, giving the vortex the same symmetry 
as the gap (Figure ||). Because the nodes are 90° apart, the superfluid velocity 
can never be perpendicular to all nodes, and there are no places where N(Ep) 
becomes zero. Both the density of states, and in a clean system the anisotropy 
in the density of states, extend throughout the entire vortex. 



B. Parks et al. p3[ have pointed out that their measurements of vortex 
dynamics in YBCO thin films are best explained by the presence of the 
extra quasiparticles expected for a vortex with an anisotropic gap. Recent 
theoretical work suggests that the structure of the vortex core in a d-wave 



superconductor is both interesting and complicated || gj, |25| . We have 
deliberately omitted the core region from Figure since we are not aware 
of any model for the core which suggests that our measurements are probing 
the core states. 



5.4 Predicted Parallel Field Dependence 

Because N(Ep, H) depends on the orientation of the current with respect to 
the nodes, improved measurements of the full angular dependence ^o tlf {H)T 
should be sensitive to the positions of the nodes. Define the orientation 
of the applied field H as (9,ip), where 9 is the polar angle measured with 
respect to the c-axis and ip is the azimuthal angle measured with respect to 
the (110) crystalline axes. For tetragonal d x 2_ y 2 symmetry with antinodes 
along the Cu-0 bonds, ip = corresponds to a field parallel to a node (Figure 
^j. In this orientation, (9 = 7r/2,yj = 0), the currents in the ab-plane are 
flowing parallel to a single node. The density of states at the Fermi level 
is y/2 lar ger at (9 = ti/2, ip = vr/4), where the currents flow parallel to 
the antinodes (Figure ^): i n this orientation, the density of states picks up 
contributions from all of the nodes. The complete in-plane angular magnetic 
field dependence of the density of states is given by 



N{0,H,9 = n/2,ip) = — (\sin<p\ + \ cos <p\)N F 



(5) 



2 M ' ' "" '^V 

as shown by Volovik ||26|| . Using 7„ = 20 mJ/molK 2 and taking H/H c2 \\ = 
0.005 as an experimentally accessible field, the field-dependent linear term 
7||(if) will vary from 0.7 to 1.0 mJ/molK 2 as a function of the azimuthal 
angle (p. This predicted angular dependence of the linear term in the specific 
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heat may be within achievable experimental resolution, although preliminary 
measurements at 8 Tesla on a heavily twinned sample with the field parallel 
to the twins and at a 45 deg angle to the twins did not detect any angular 
variation in 711 (H). 

6 Summary 

The residual density of states at the Fermi level in the superconducting state 
of single-crystal YBa 2 Cu307_5, N(Ep,H = 0)/N n , is determined by the 
zero-field linear term ^(H = 0)T in the specific heat. N(Ep,H = 0)/N n 
is qualitatively larger than would be expected from fitting the temperature 
dependence of the penetration depth [0 to expressions for lines of nodes with 
scattering |J. In two twinned single crystals, ^{H = 0)/j n ~ 0.15, while in 
two untwinned single crystals, ^(H = 0)/j n ~ 0.10. 

The specific heat also includes a j± (H)T term, which obeys 7_i_(i?) ~ 
j n \J H / H c2 as predicted for superconductors with lines of nodes in the gap 
function ||. This nonlinear field dependence of the density of states sug- 
gests quasiparticle excitations outside the vortex core, and appears to be 
independent of twinning, unlike the zero-field linear term. The in-plane an- 
gular magnetic field dependence of the specific heat, which is smaller than 
the perpendicular magnetic field dependence, is predicted to be sensitive to 
the locations of the nodes. 

We thank M. R. Beasley, C. M. Varma, and especially G. Volovik for 
many useful discussions, and R. A. Fisher, N. E. Phillips, and J. Gordon for 
sharing their analysis prior to publication. 
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Figure 1: Total specific heat at and 8 Tesla for sample Tl 



Figure 2: Nonphonon specific heat for samples Tl and Ul, with the global 
fit described in the text. 



Figure 3: Coefficient of the linear term for crystals Tl and Ul, determined 
by the global fit described in the text. 



Figure 4: Coefficient of the linear term for crystal Tl, determined by a) 
independent fits in which all parameters might be field- dependent, and b) a 
global fit with the H = data set excluded 



Figure 5: Schematic of the local density of states N(E F )/N n throughout a 
d-wave vortex. 



Figure 6: The azimuthal field angle (p relative to a tetragonal d-wave order 
parameter, and sketch of the current flows associated with an in-plane vortex 
for ip — and <p — ir/4. 
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